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^ ■ Abstract 

^ • Given a Boolean function f on n variables, a Disjoint Sum-of-Products (DSOP) of / is a set of 

products (ANDs) of subsets of literals whose sum (OR) equals /, such that no two products cover the 



same minterm of /. DSOP forms are a special instance of partial DSOPs, i.e. the general case where 
I a subset of minterms must be covered exactly once and the other minterms (typically corresponding to 

don't care conditions of /) can be covered any number of times. We discuss finding DSOPs and partial 
DSOP with a minimal number of products, a problem theoretically connected with various properties 
of Boolean functions and practically relevant in the synthesis of digital circuits. Finding an absolute 
minimum is hard, in fact we prove that the problem of absolute minimization of partial DSOPs is NP- 
hard. Therefore it is crucial to devise a polynomial time heuristic that compares favorably with the 
I known minimization tools. To this end we develop a further piece of theory starting from the definition 

of the weight of a product p as a functions of the number of fragments induced on other cubes by the 
selection of p, and show how product weights can be exploited for building a class of minimization 
heuristics for DSOP and partial DSOP synthesis. A set of experiments conducted on major benchmark 
functions show that our method, with a family of variants, always generates better results than the ones 

■ of previous heuristics, including the method based on a BDD representation of /. 

m ■ 
in 

■ 1 Introduction 

O 

CN [ Given a Boolean function f on n variables x\,X2,...,Xn in "S", a Disjoint Sum-of-Products (DSOP) of / is 

a set of products (ANDs) of subsets of literals whose sum (OR) equals /, such that no two products cover 
the same minterm of /. As each product is the mathematical expression for a cube in S", a DSOP also 
^ I represents a set of non intersecting cubes occupying the points of "S" in which / = 1. In fact we shall 

^ ■ indifferently refer to products or cubes, and apply algebraic or set operations to them. We are interested in 

- - - finding a DSOP with a minimal number of products. 

Besides its theoretical interest, DSOP minimization is relevant in the area of digital circuits for deter- 
mining various properties of Boolean functions and for the synthesis of asynchronous circuits, as discussed 
for example in ||4lll0j|TT][l2j[T6]. DSOPs are indeed used as a starting point for the synthesis of Exclusive- 
Or-Sum-Of-Products (ESOP) forms, and for calculating the spectra of Boolean functions. 

DSOP forms can be seen as a special case of partial DSOPs where a subset of minterms of a Boolean 
function must be covered exactly once, while other minterms can be covered more than once or not be 
covered at all. In particular this is the case where the points in the on set of a function are covered exactly 
once, while the points in the don't care set can be covered any number of times [ 10|. 

For speeding an otherwise exceedingly cumbersome process an absolute minimum in general is not 
sought for, rather heuristic strategies for cube selection have been proposed, working on explicit product 
expressions ||2j|5]|T5l, or on a BDD representation of / [3] IS. 
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After discussing the complexity of DSOP and partial DSOP absolute minimization we propose a class of 
heuristic algorithms based on the new concept of "cube weight", and show that our results compare favorably 
with the ones of the other known heuristics. The starting set of cubes is the one of a sum of product (SOP) 
found with standard heuristics. The SOP cubes may be eventually fragmented into non overlapping sub- 
cubes, giving rise to a largely unpredictable DSOP solution. The process may exhibit an exponential blow 
up in the number of fragments even dealing with theoretically minimal solutions, as for a function presented 
in fB| where |SOP| = n/2 and |DSOP| = 2"/^ - 1 (|SOP| and |DSOP| denote the number of terms in the 
SOP and DSOP expression, respectively). 

Another new characteristic of our heuristic is the idea of recomputing a SOP on the residual function 
at different possible stages of the disjoint minimization process, as a tiade-off between quality of the result 
and computational time. We have observed experimentally that this strategy is crucial for obtaining compact 
DSOP forms. For ease of presentation we start with DSOP synthesis and then extend the heuristics to the 
more general case of partial DSOP. 

The paper is organized as follows. In the next Section |2] we discuss the complexity of absolute min- 
imization of DSOP and partial DSOP forms proving that, for the latter, i.e. for the most general forms, 
the problem is NP-hard. In Section |3] we define the weight of a product as a function of the number of 
fragments possibly induced on other cubes by the selection of p. In Section |4] we show how this weight can 
be exploited for building a class of minimization heuristics. Section |5] extends our strategy to partial DSOP 
synthesis. In Section [6] we present and discuss the computational results obtained by applying the proposed 
heuristic to the standard ESPRESSO benchmark suite I.18J . and comparing these results with other published 
data. The paper is concluded in Section |7] 

2 The complexity of DSOP minimization 

As it may be expected absolute DSOP minimization is a hard problem and absolute partial DSOP minimiza- 
tion may be at least as hard. Let us first recall some classical definitions. In a Boolean space {0, 1 }" described 
by n variables x\,X2, x„, a completely specified Boolean function is a function / : {0, 1}" — > {0, 1}, while 
Boolean a function / is partial if / : {0, 1}" {0, 1, — }. With usual terminology, a literal yi is a variable 
Xi in direct or complemented form, and products are ANDs of literals. A product p is an implicant of the 
Boolean function / if Va; G {0, 1}", {p{x) = 1) ^ ifi^) = 1)- An implicant p of a. function / is a prime 
implicant if p cannot be implied by a more general (i.e., with fewer literals) implicant of /. 

Unlike SOPs, a DSOP composed of prime implicants only may not exist, as can be immediately seen 
considering a function with only three points in the on set, one adjacent to the other. Furthermore, DSOPs 
of prime implicants may exist but none of them may be minimal. For example the minimal DSOP cover 
of six implicants shown in Figure 1 contains the non prime implicant x\Xd,xiX(pci displayed in the sub-map 
xsX(^-i = 001, which is covered by the prime implicant xiX4X5X(, spanning across the sub-maps X5X(,xj = 
000 and xsx^xy = 001. The reader may discover that there is one DSOP cover composed of seven prime 
implicants but not less (actually we could not construct an example with less than seven variables). 

The above considerations show that, unlike in the SOP case, in DSOP minimization non prime impli- 
cants must also be considered. Theoretically this is not a major drawback as the generation of all implicants 
requires polynomial time in the size of the input (truth table of the function). The problem arises in the impli- 
cant selection phase where, as in the SOP case, a brute force enumerative selection requires exponential time 
in the worst case. It has been shown that SOP absolute minimization is as complex as set covering iTTlfTTl. 
Similarly DSOP absolute minimization can be compared to the set partitioning (or minimal exact cover) 
problemQ It is immediate that minimal exact cover is at least as hard as absolute DSOP minimization (solv- 

'The minimal exact cover problem is as follows: given a family of subsets 5 of a set U and a positive integer k, is there a subset 
family T CS such that the subsets in T are k in number, are disjoint, and their union is the entire set t/? The minimal exact cover is 
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Figure 1: A minimal DSOP in seven variables composed of six cubes, one of which is not prime (the 
Karnaugh maps for X5X6X7 =011, 101, 110, 111 do not contain I's and are not shown). A corresponding 
DSOP composed of prime cubes only includes at least seven of them. 



ing the former problem efficiently would imply solving also the latter). Here we are not proving the reverse 
condition, rather we focalize on the most general problem of partial DSOP absolute minimization and prove 
that, in this version, the problem is NP-hard. More precisely, we prove that the decision version of partial 
DSOP minimization is NP-complete. Let us formally define the problem. 

MIN Partial DSOP 

Input: A partial Boolean function / : {0, 1}" — ?• {0, 1, — }, specified by its on, off, and don't 
care set, and a positive integer k. 

Question: Is there a partial DSOP, i.e., a sum of products covering exactly once the points of 
the on set, and any number of times the points in the don't care set of /, with at most k products? 

This problem is in NP because given a candidate partial DSOP with at most k terms, one can determine 
whether it is a covering of / satisfying the given requirements in time polynomial in the size of the input 
instance. In fact this simply requires evaluating the partial DSOP at all of the points in the on set of / and 
checking that one and only one of its products takes the value 1 . 

To prove the NP-completeness of MIN Partial DSOP, we adapt to our problem the theory and the 
proofs developed in [1], where the authors proved that the decision version of finding the smallest SOP form 
consistent with a truth table is NP-complete, reducing from 3-Partite Set Cover, instead of Circuit 
SAT as done in fTTIl . Moreover, they pointed out that the reduction would also work for 3D MATCHING, 
which is precisely the NP-complete problem that we will reduce to MiN PARTIAL DSOP. 

3D Matching 

NP-hard since it can be easily reduced by the "exact cover problem" introduced by Karp in 1972 f9l setting k to the cardinality of 
U. 
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Input: A positive integer n, a partition IT of the set {1,2, ... ,n} into three sets of equal size, 
and a collection S of subsets of {1,2, .. . ,n}, where every subset contains exactly one element 
from each of the set of n. 

Question: Is there a subcollection C C 5 of size n/3 whose union is {1,2, .. . ,n}7 

Note that such a subcollection C would provide an exact cover, as it covers each element of the set { 1 , 2, . . . , «} 
exactly once. 

In the next theorem, we give the reduction from 3D MATCHING to MiN PARTIAL DSOP. It is basically 
the same reduction given in HI for SOP minimization, however here we show how it works even for disjoint 
SOP minimization. Given m, v G {0, 1}", we will write m < v if m,- < v,- for all / G {1,2, . . . ,«}. 

Theorem 1 MiN partial DSOP is NP -complete. 

Proof. We have already noticed that MiN PARTIAL DSOP belongs to NP. Thus, we are left to show that 
it is NP-hard. To this aim we show how to transform an input instance of 3D MATCHING into an instance 
of MiN PARTIAL DSOP in polynomial time. The instance defines an incompletely specified function / 
depending on O(logn) variables, that can be covered by a partial DSOP with n/3 products if and only if 
there is a subcollection C C J of size n/3 whose union is {1,2, ...,«}. 

Let (n,n,5) be an input instance of 3D MATCHING. We first define two sets of vectors, V and W, 
that we will use to define an instance of MiN PARTIAL DSOP. Let q be the smallest even integer such that 
( %) ^ Observe that q = C?(log«). We assign a unique q-bit vector b{i) with exactly q/2 I's to each 



/ G {1,2, . . . ,«}. Let n(/) G {1,2,3} be the index of the block of the partition 11 that contains /. Let t = 3q. 
The vectors in V and W can be divided into 3 blocks, each of size q. We can now define the vectors in 

V = {v(') \\<i<n}andW = {w^^' | A G 5}: 

• each v^'' G V , / G { 1 , 2, . . . , «}, is equal to b{i) on block n(/), and is in the other two blocks; 

• each w^'^' G IV, A G 5, is the bitwise OR of all v^'' G V such that / G A. 

These two sets can be generated in time n'^^^\ 
Observe that this choice guarantees that: 



The forward implication is obvious. To see that the backward implication holds, let A G 5 and / g{1,2,...,«}, 
and assume that v^'^ < w^^\ This implies that A contains one element j that belongs to the same block n(/) 
of /, where v^'^ is not 0, i.e., n(/) = 11(7). Thus, since v^'^ < w^^\ we must have b{i) < b{j), which in turn 
implies / = j, and therefore / G A. 

We now construct an incompletely specified function / on the domain {0, l}^ as follows: 

• f{x) = 1 ifxG V. 

• f{x) = — ifx^V and .jc < w for some w G W. 

• f{x) = 0, otherwise. 

For u G {0, 1}', let D{u) = {w \ w < u} and let x(m) denote the product Yli:ui=o^i- Observe that t(m) is 
the characteristic function of the set D{u). Consider the set D{W) = [J^^^,D{x). Property (1) implies that 
V C D{W) and that f{x) = - iff ;c G D{W) \ V. 

To complete the proof we must show that S contains a cover C of size n/3 if and only if there is a partial 
DSOP for /, with n/3 products. Suppose that S contains a cover C of size n/3, and consider the set of 
products {x(w('')) I C G C}- It is immediate to verify that the sum of these products covers the function 





n} 




(1) 
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/. Indeed, for all / G {1,2, . . . ,n}, i belongs to one of the sets in C, say C' , and the vector v(') in the on 
set of / is then covered by the corresponding product x(w('' ^) (recall that by consti^uction v^'^ < w'^^\ thus 
v(') € D(w^^^)). Now, we have to prove that these products define a partial DSOP for /, i.e., we must show 
that the corresponding cubes are either disjoint or intersect only on the don't cares of /. 

First of all recall that each vector w^^\ C ^ C, can be divided into three blocks, each equal to one of 
the vectors b{i). For instance, if C = {i,j,k\, with n(/) = 1, n(j) = 2, and n(^) = 3, then w^''^ is given 
by the concatenation of b{i), b{j), and b{k). The related product 'z{w^^^) can then be divided into three 
subterms of q/2 literals, containing the complemented variables corresponding to the O's in b{i), b{j), and 
b{k). Moreover, since C is a disjoint cover of {1,2, . . . ,n}, each v^'^ belongs to one and only one of the sets 
in C, that is only one of the vectors w^'^^ has a block equal to b{i), for all / G {1,2, . . . ,«}. This implies that 
all subterms of the set of products {x(w^'^^) \ C £ C} are different. 

Given any pair of products x(w(*-^') and x(w(^'), with C,D G C, consider the intersection of the corre- 
sponding cubes. The characteristic function of the intersection is simply the product (AND) between x(w(*'') 
and x(h''^'). Since all subterms of z{w^^^) and x(w(^)) are different, the product x(>v'^''^) •x(w^^^) contains 
three subterms, each of at least q/2 + 1 complemented variables. Thus, it can cover only don't cares of /, 
since any vector v''' G V has one block with only q/2 O's. 

Now, suppose that (j) is a partial DSOP for /, with n/3 products. For each product p let u{p) be the 
maximal vector satisfying p. Note that G {I,—}, thus u{p) £D{W) and there must be a set S(p) g5 
such that u{p) < w^^^p^\ We then show that the collection C = {S{p) | p G (|)} is a cover of {1,2, . . . ,«}. Let 
7 G {1,2, . . . Since f(y^-'^) = 1, exactly one of the product in (|), say p^^\ must cover v'^'. This implies 
v^-/) < u{p^}^). Thus v(^) < w^^^P^'^^K which by property (1) implies ; G S{p^j'^). ■ 

The exponential nature of partial DSOP minimization justifies the search for heuristic solutions. This 
will be done after a theoretical discussion on how cubes gets fragmented due to their intersections, contained 
in the next section. This will lead to a heuristic strategy whose complexity is polynomial in the size of the 
output, i.e., in the number of products of the computed DSOP form. 

3 The Weight of a Cube 

A product q = yi^yh...yi|., \ <k <n, represents a cube of dimension d{q) = n — k, i.e., a cube of 2"^*^ points 
in {0, 1}". The intersection p = pinp2of two cubes p\ = yi^...yi|^^ , p2 = jj, ...j,^^^ is obviously obtained as 
the AND of the two corresponding products. The intersection p is empty if and only if there is a literal in 
pi that appears complemented in p2, and vice-versa. Otherwise p is a cube of dimension d{p) = r, with 
r = n — {k\+k2 — c), and c is the number of common literals in p\ and p2. 

Take p\, p2 as above, and let p\, p2 partially overlap. The set of points of P2\pi can be covered in 
different ways by a set of at least ki —c disjoint cubes of dimensions r,r+l,...,?i — ^2 — 1- For n = 6, letting 
ki = 5, k2 = 3, c = 2 we have r = and d{p\) = 1, d{p2) = 3, i.e., the intersection contains 1 point, and the 
two cubes contain 2 and 8 points, respectively. Therefore, P2\pi contains 7 points and can be covered with 
5 — 2 = 3 cubes of dimensions 0, 1, 2. For an other example, consider cubes A and B in Figure |2ta). The set 
A\B contains the minterms 0000, 0001, and 0100. The disjoint covers for these points are X1X2X3 +x\X2X^M 
and X1X3X4 +xiX2^3^4, both containing two cubes. 

Now, if p\ is selected into a DSOP, p2 must be discarded and the points of p2 \pi must be covered 
with at least k\ — c disjoint cubes instead of one (the single P2). Then k\— c — \ is the number of exti^a 
cubes required by the DSOP. If the function / can be represented by a SOP containing only p\ and p2, the 
selection of p\ into a DSOP requires a total of A;i — c + 1 cubes. In particular if ^1 — c = 1 the intersection p 
covers exactly one half of the points of p2 and P2\pi is also a cube. Clearly the general situation will not 
be that simple as the starting SOP for /, to be transformed into a minimal DSOP, will consist of a collection 
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Figure 2: (a) A minimal SOP of four cubes of dimension 2 in "3^, with weights w(A) = 1, w{B) = 2, 
w(C) = 0, w(D) = 1. (b) A corresponding DSOP. 



of cubes overlapping in groups. Still we define a weight for each cube pi equal to the minimum number of 
extra cubes that the selection of pi would induce in all the cubes intersecting p,. Formally, let a SOP for / 
consist of partially overlapping products pi,p2,...,Ps- We pose: 

Definition 1 Let a product pi of k literals intersect the products pi^,...,pii, such that pi and pi have cj 
common literals. Then w{pi/ pi-) = k — cj — I is the weight of pi relative to pj., and w{pi) = L;=i ^{Pi/ Pij) 
is the weight of pi. If pi does not intersect any other product, set w{pi) = —I. 

Thus, when pt intersects , the weight of pi relative to is the minimum number of additional 
products that we would have in the cover keeping pi and covering pi/ pi. with non-overlapping products. 

As an example, consider the function / of four variables, represented in Figure 0a). A minimal SOP 
of / contains four cubes A = I1X3, B = X2X4, C = x\X2, D = x^xt,, all of dimension two. The weights 
are computed as follows. For A: w{A/B) = 1 (in fact, selecting A in a DSOP would require to cover the 
remaining three points of B with at least two disjoint cubes); w(A/C) = (the residual two points of C can 
be covered with one cube); then w(A) = 1. For S: w{B/A) = 1; w{B/C) = 0; w{B/D) = 1; then w{B) = 2. 
For C: w(C/A) = 0; w{C/B) = 0; then w(C) = 0. For D: w{D/B) = 1; then w{D) = 1. As we shall explain 
in the next section, we start the construction of a DSOP by selecting the cubes with low weight and high 
dimension, breaking on the fly the ones that intersect a selected cube. In the present example, start by 
selecting C and reduce A and B to two subcubes Ai, Bi of two points each. Then select D and further reduce 
Bi to B2 of one point. Then select A\ and B2, as shown in the DSOP of Figure^b). During the process the 
weights are updated as explained below. 

4 DSOP synthesis algorithms 

Let us consider an incompletely defined Boolean function / : {0, 1}" — > {0, 1, — } represented with a set of 
cubes C = {Con,Cdc), where Con covers the on set of /, i.e., the points v in {0, 1}" such that /(v) = 1, and 
Cdc covers the don't care set of /, i.e., the points v in {0, 1}" such that /(v) = — . 

The new heuristic for DSOP construction uses four basic procedures working on an explicit represen- 
tation of cubes. The first procedure BUILD-SOP(C,P) works on a set C of cubes covering an arbitrary 
function as above, to build a minimal (or quasi minimal) SOP P for that function. Note that, during the 
process, BUILD-SOP may be called on different sets C emerging in the computation. As a limit the cubes 
of C may be minterms, i.e., cubes of dimension 0. The second procedure WEIGHT(P) builds the weights 
for the cubes of a set P. 
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algorithm DSOP(C,D) 

INPUT: A set of cubes C covering a function / 
OUTPUT: A set of disjoint cubes D covering / 

D = % 

while (C ^ %) 

BUILD-SOP(C,P) 

A = {dcP\\lccP\{d}: dC\c = %} 

D = DUA 

P = P\A 

WEIGHT(P) 

SORT(P) 

B = 

while (P ^ 0) 

let p be the first element of P 

P = P\{p} 

D^DU{p} 

forall qeP: pCiq^Q) 

P^P\{q} 

BREAK(^,p,0 

OPT(^,e,P,B) 
forallreB: pHr^Q 

B = B\{r} 

BREAK(r,p,Q) 

B = BUQ 

C^B 

Figure 3: The general algorithm for DSOP synthesis. 

The third procedure SORTCP) sorts a set P of weighted cubes. This procedure comes in two versions: i) 
the cubes are ordered for decreasing dimension and, if the dimension is the same, for increasing weight; ii) 
the cubes are ordered for increasing weight and, if the weight is the same, for decreasing dimension. If two 
or more cubes have same weight and same dimension, their order is chosen arbitrarily. The two versions of 
SORT give rise to two different alternatives of the overall algorithm. 

The fourth procedure BREAK{q,p, Q) works on the set difference q\p between two cubes, to build an 
arbitrary minimal set Q of disjoint cubes covering q\p. Note that this operation is easy since q\p can be 
obtained as q\{pr\q), where the latter is the set difference between two cubes, i.e., q and pHq, in turn a 
cube because is the intersection of two cubes. 

In practice, for BUILD-SOP one can use any minimization procedure (in our experiments we have used 
procedure ESPRESSO-NON-EXACT of the ESPRESSO suite QH). Procedures WEIGHT and SORT (both 
versions) are obvious. Procedure BREAK is the one suggested in f8| and 1 14] as DISJOINT-SHARP. 

In the overall process we consider four sets of cubes C,P,B,D. At the beginning C contains the cubes 
defining /, while P,B,D are empty. During the process C contains the cubes defining the part of / still to 
be covered with a DSOP; P contains the cubes of a SOP under processing; B temporarily contains cubes 
produced by BREAK as fragmentation of cubes of P; and D contains the cubes already assigned to the 
DSOP solution and, at the end, the solution itself. 

The algorithms of our family share the structure shown in Figure |3] (its behaviour on incompletely 
specified functions is discussed at the end of this section). As long as / has not been completely covered 
with disjoint cubes, i.e., there are still cubes in the set C, a minimal (or quasi-minimal) SOP P for the part 
of / still to be covered is computed by the procedure BUILD-SOP. All cubes that do not intersect any other 
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cube in P are removed from P and inserted in the DSOP D under construction; the remaining cubes are 
weighted and sorted. Then, the first cube p is extracted from P and inserted in the solution D. Each cube 
q G P that intersects p is removed from P, and a SOP Q for the set difference q\p is computed by the 
procedure BREAK. 

During this phase an optional optimization procedure OPT is called to decide how to handle the frag- 
ments in Q; depending on this optimization phase, different variants of the heuristic can be defined. Note 
that, since the points of p cannot be covered by any other cube, all fragments r already inserted in B must 
be tested for intersection with p and, if necessary, replaced with the SOP computed by BREAK for the set 
difference r\p. When P becomes empty, the fragments in B are moved to the set C and the algorithm iter- 
atively builds a new SOP P covering the points that are not yet covered by the DSOP D under construction. 
The iterations terminate when C becomes empty. 

We have designed and tested five variants of our heuristic based on five different versions of the opti- 
mization procedure OPT, with different degrees of sophistication. The first variant, DSOP-1, is the simplest, 
and computationally fastest, as OPT simply inserts the cubes of Q into the set of fragments B: 

in DSOP-1: 

procedure OPT(q,Q,P,B) 
B = BUQ 

Example 1 For an example, Figure W[b) shows a DSOP form for the SOP form of Figure ^a), com- 
puted by algorithm DSOP-1. At the beginning D = % and P = {xiX2,xiX3,xiX3,X2X4}, sorted for de- 
creasing dimensions of cubes and then for increasing weights (we recall that, w(^\X2) = 0, w{x\x^) = 1, 
w{^{x^) = 1, and w{x2X4) = 2). The first cube considered is p = x\X2, which is removed from P and in- 
serted in D. Its intersecting cubes, x\Xt, and X2X4, are then broken generating the residuals cubes xixjx^ 
and X1X2X4, respectively, which are inserted in B, while X1X3 and X2X4 are removed from P. The last cube 
in P to be considered is then x\x^ that is inserted directly in D, since there are not any other remain- 
ing cubes in P. Its intersecting cube X1X2X4 in B is then reduced to xyx^x^Xd,. The second while (P ^ 0) 
iteration starts with P = {xiX2X^ ,X[X2XjX4} and D = {x[X2,xiXt,}, and terminates with the final DSOP 

D = {xiX2,XiX3,XiX2X3,XiX2X3X4}. 

In the second variant, DSOP-2, after a cube p has been selected and moved to D, each cube q intersecting 
p is, as before, fragmented and moved to B. In addition the optimization procedure updates the weight of 
all cubes r e P that intersect q, and then sorts the cubes in P again: 

in DSOP-2: 

procedure OPT{q,Q,P,B) 
B = BUQ 

/ = {r €P I qnry^d)} 

WEIGHT(/) 

SORT(P) 

A disadvantage of both versions is that whenever a cube p is moved from P to D, all cubes q intersecting 
p are fragmented and removed from the set P. Hence, the fragments, even the big ones, are "out of the 
game" and cannot participate in the construction of the DSOP D until P becomes empty and a new SOP 
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covering all fragments in the set B is computed. Consequentially, small cubes in P could be selected first, 
possibly damaging the quality of the final result, i.e., the size of the final DSOP. 

To partially avoid this disadvantage, we have implemented a third version of the heuristic, DSOP-3, in 
which whenever a cube p G P is moved to D, each cube q intersecting p is, as before, fragmented and moved 
to B, and, in addition, all cubes r G P intersecting q are moved to B as well: 



in DSOP-3: 

procedure OVTiq,Q,P,B) 

/= {r GP I qOr^^ 
B = BUQUI 
P = P\I 



In this way, the cubes of P intersecting the fragments already in B cannot be selected, while is avoided 
the possible fragmentation of big cubes in B. Moreover, we leave open the possibility of selecting these big 
cubes in the next iterations of the algorithm. This version of the heuristic is computationally more expensive, 
since in the internal while loop less cubes can be selected (P empties faster), and procedure BUILD-SOP 
must be executed more frequently. 

The fourth version of the heuristic, DSOP-4, checks whether the set Q contains only one fragment, i.e., 
^ \ is a cube. In this case, this only fragment is put back in P. The cubes left in P are then weighted and 
sorted again: 

in DSOP-4: 

procedure OPT(q,Q,P,B) 

if(iei = i) 

P = PUQ 
else 

B = BUQ 
WEIGHT(P) 
SORT(P) 



Finally, in the last version of the heuristic that we have tested, DSOP-5, the biggest fragment in the set 
Q is always put back in P. The cubes left in P are then weighted and sorted again. In this way, big fragments 
remain part of the game in the present iteration of the algorithm: 

in DSOP-5: 

procedure OPT(q,Q,P,B) 

let b be the biggest cube in Q 

P = PU{b} 

B = BUQ\{b} 

WEIGHT(P) 

SORT(P) 
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The performances of these five procedures are discussed in Section[6] We have observed experimentally 
that more sophisticated optimization procedures do not always provide better quality results. Experimental 
results have also outlined how the BUILD-SOP procedure, i.e., re-synthesizing the remaining cubes, seems 
to be crucial for obtaining compact DSOPs. 

Let us now briefly consider the case of the DSOP synthesis of incompletely specified Boolean functions. 
Our heuristic does not consider explicitly the presence of don't cares; indeed, the first call of the BUILD- 
SOP procedure produces a SOP P covering the whole on set of / and a subset of its don't care set. Then, 
the algorithm works on the SOP P, treating all points covered by its cubes as if they belonged to the on set 
of /, i.e., there is no distinction between points originally in the on set of / and points originally in the don't 
care set. In particular, the successive calls of BUILD-SOP on the part of / still to be covered with a DSOP, 
treat the function as if it were completely specified. Of course, each cube in the SOP P computed by the 
first call of BUILD-SOP covers at least one point in the on set of /, as cubes covering only points in the 
don't care set are discarded by the SOP minimization algorithm. However, the final disjoint cover D for / 
could contain cubes covering only points originally in the don't care set. In fact, cubes in D are either entire 
cubes of the starting SOP P, or sub-cubes of cubes in P (besides new cubes and sub-cubes originated by the 
successive calls of BUILD-SOP) and some sub-cubes (or new cubes) could only cover don't care points. 

From the above all versions of our heuristic could be improved checking whether a cube p contains only 
points in the don't care set of the function /, before adding it to the DSOP solution D under construction. 
Unfortunately, such a check can be computationally expensive, and for this reason we have not added it as 
a "default" procedure in our algorithm. In fact the check is left as an option. Experiments conducted on a 
set of incompletely specified functions show some improvements on the final form induced by the check at 
a considerable increase of computing time, see next Section |6] 

5 Partial DSOP synthesis 

As already mentioned the problem of DSOP minimization naturally generalizes to covering partial DSOPs 
where some minterms (e.g. the ones in the on set of the function) are covered exactly once while other 
minterms (e.g. the ones in the don't care set) can be covered any number of times lITOl . In this section we 
present a general heuristic to efficiently compute a partial DSOP cover. 

The heuristic makes use of two sums of products as input. The first SOP, sopD, contains all points of 
the on and don't care set of the function / that must be covered only once (DSOP part), while the second 
SOP, sopS, contains all the points of / that can be covered more than once (SOP part). These two SOPs 
are disjoint. The output of the heuristic is a cover of the overall function /, represented by the union of the 
two SOPs sopD and sopS that respects the specifications. Note that when sopD is empty the problem is a 
classical SOP minimization, while when sopS is empty the problem is a classical DSOP minimization. 

The algorithm uses four basic procedures as for the DSOP synthesis of SectionHl In particular BUILD- 
SOP, WEIGHT, and SORT are the same. 

The fourth procedure FARTlAL-BREAKiq, p,sopD,sopS,Q,R) works on the set difference q\p be- 
tween two cubes, to build an arbitrary minimal set Q of disjoint cubes covering q\p, if qCi p is not entirely 
contained in sopS. If q Hp is contained in sopS, the cube q is not broken and we can keep it in the set P 
which contains the cubes to be considered in the current iteration. In this case we then set 2 = 0. Moreover, 
the procedure PARTIAL-BREAK builds a set R containing points of q\p that can be covered more than 
once and can therefore be added as don't cares to C. In this way, these points, that have been already cov- 
ered, could be used again in the minimization phase to get a smaller cover. This procedure, different from 
the one used for DSOP synthesis, is presented in Figured 

The overall minimization heuristic is presented in Figure |5] As for the DSOP synthesis, the heuristic 
makes use of four sets of cubes C,P,B,D. At the beginning C = sopDU sopS contains the cubes defining 
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algorithm VARTlAL-BT(EAKiq,p,sopD,sopS, Q,R) 
INPUT: The chosen cube p, the cube q that can be broken, 
the two SOPs sopD and sopS whose union represents / 
OUTPUT: A minimal set Q of disjoint cubes covering q\p and 
a set R of the points of q\p that can be covered more than once 

R = <d 

Pi=qnp 

if (Pi Q sopS)// all points of pi can be covered more than once 
Q = <b 

else if (pi C so pD) // all points of pi must be covered once 

Q = DISJOINT.SHARP(^,p,) 
else // Pi intersects both sopD and sopS 

Q = DISJOINT.SHARP(^,/70 

R = PiCi sopS 

Figure 4: The procedure PARTIAL-BREAK to be used in partial DSOP synthesis. 

algorithm FARTlAL-T>SOPisopD,sopS,D) 

INPUT: Two disjoint SOPs describing the points of / that 

must be covered only once (sopD) and the points of / that 

can be covered more than once (sopS) 

OUTPUT: A partial DSOP D for the function / 

Con = SOpDon UsopSgn 

Cdc = sopDdc U sopSdc 
whUe {Con ^ 0) 

BUILD-SOP(C,P) 

A = {d eP\'ic&P\{d}: d^c = <b] 

£) = DUA 

P = P\A 

WE1GHT(P) 

SORT(f) 

5 = 

while {P ^ 0) 

let p be the first element of P 
P = P\{p} 
D = DU{p} 
foraHq&P : pCiq^d) 

FARTlAL-BREAK(q,p,sopD,sopS, Q,R) 

it(Q^&)P = P\{q} 

OVT(q,Q,P,B) 

Cdc = Cdc^R 
forallrGB: pdr^d 

PARTlAL-BREAK(r, p, sopD, sopS, Q,R) 

if(e/0)fi = B\W 

B = BUQ 

Cdc = Cdc^R 

Con =B 

Figure 5: Algorithm for partial DSOP synthesis. 
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Figure 6: (a) sopS (cubes with solid lines) and sopD (cubes with dotted lines), (b) A corresponding partial 
DSOP. 



/ while P,B,D are empty. During the processing C contains the cubes defining the part of / still to be 
covered with a partial DSOP. P contains the cubes of a SOP under processing. B temporarily contains cubes 
produced by BREAK as fragmentation of cubes of P. D contains the cubes already assigned to the partial 
DSOP solution and, at the end, the solution itself. OPT{q,Q,P,B) is an optional optimization procedure 
to decide how to handle the fragments produced by the procedure BREAK. As before, depending on this 
optimization phase, different variants of the heuristic can be defined. 

Example 2 Consider the function shown in Figure^a). Suppose that sopD = {x[X2X3,xiX2XjX4} (cubes 
with dotted lines in the figure) and sopS = {xiXj,xiX3} (cubes with solid lines). A partial DSOP for f is 
shown in Figure ^b). This expression is obtained with the partial DSOP algorithm as described in the 
following. Let OPT(q,Q,P,B) be the simple command B = BUQ (as in the DSOP-1 procedure). At the 
beginning D = % and, after the SOP minimization phase, P = {x\X2,x\Xt,,x\Xt,,X2X4}, sorted for decreasing 
dimensions of cubes and then for increasing weights ( note that we have the same initial P of Example [7]). 
The first cube p = x\X2 is removed from P and inserted in D. Its intersecting cubes are Xix^ and X2X4. In the 
procedure PARTIAL-BREAK, the intersection between X2X4 and x\X2 is pi = X1X2X4. Note that pi intersects 
both sopD and sopS, thus Q = {X1X2X4} and R = {^1X2X3X4} (i.e., X1X2X4 and X1X2X3X4 will be inserted in 
B and in the don't care set ofC, respectively). Moreover, we compute the intersection pi between Xixj and 
X1X2, obtaining X1X2X3 which is entirely contained in sopS. Thus, in this case Q = R = %, then x^x^ is not 
broken and it is not removed from P. Similar operations are performed on X1X3, and on X1X2X4 contained 
in B. The second while (P 7^ 0) iteration, which starts with the P = {X1X2X4} and D = {xiX2,xiX3,xiX3}, 
terminates with the partial DSOP shown in Figure^ b). 

6 Experimental Results 

In this section we present and discuss the results obtained with the heuristics presented above to the standard 
ESPRESSO benchmark suite [18|. All experiments were performed on a 1.8 GHz PowerPC with 1 GB of 
RAM. 

6.1 DSOP synthesis 

We have considered the five different variants of the heuristic described in Section IH denoted as DSOP-1, 
DSOP-2, DSOP-3, DSOP-4, DSOP-5. For each variant, we have run both versions of the procedure SORT, 
to estimate the practical effectiveness of each version. Namely we have ordered the cubes for decreasing 
dimension and, in case of equal dimension, for increasing weight (version dimension/weight). Then we 
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0.99 
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32 


20 
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1.36 
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0.84 
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1.33 


303 


0.54 
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0.56 
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15 


7 
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811 
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14 


14 
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2.57 
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94 
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0.59 
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10 


1 
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0.43 
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0.51 
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1.11 
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0.85 
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14 


14 
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0.41 
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0.23 
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0.33 


179 


0.16 
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0.18 


table5 


17 


15 


158 


167 


0.39 
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0.36 


161 


0.38 


161 


0.24 


161 


0.25 


tial 


14 


8 


581 


943 


1.78 


937 


1.96 


874 


2.98 


1071 


2.84 


1040 


2.50 


vtxl 


27 


6 


110 


204 


0.45 


204 


0.49 


204 


0.64 


208 


0.34 


213 


0.31 


x7dn 


66 


15 


538 


796 


1.30 


784 


1.43 


812 


1.57 


813 


0.88 


864 


0.76 



Table 1: Compaiison of five different valiants of the DSOP minimization heuristic (SORT version: dimen- 
sion/weight.) The size of the best DSOP representation computed for each benchmark is in boldface. 

have ordered the cubes for increasing weight and, in case of equal weight, for decreasing dimension (version 
weight/dimension) . 

Since the benchmarks are multi-output functions and the algorithm is described for single output func- 
tion, in the experiments we have considered each output separately, but the minimization phase with ESPRESSO 
is performed in a multi-output way. Moreover, common disjoint cubes of several output are counted only 
once. 

Tables [U and |2] report a significant subset of the experiments. In particular. Table [T] reports the perfor- 
mances of the heuristics with respect to the first version of the SORT procedure, while Table |2] is relative 
to the second SORT procedure. All benchmarks in these tables are completely specified. In both tables, 
the first column reports the name of the benchmark; the following two columns give the number of inputs 
and outputs; the column labeled SOP shows the number of products in a SOP representation computed by 
ESPRESSO in the heuristic mode; finally the remaining five pairs of columns report the number of disjoint 
products in the DSOP expressions computed by our heuristics and the corresponding synthesis time. 

As Table [T] and Table |2] clearly show, the third variant of the heuristic, together with the first version 
of procedure SORT (version dimension/weight), gives the best results regarding the size of the resulting 
DSOP forms, and its running times are comparable to those of the other variants, and sometimes even lower. 

We have then tested the performances of the best variant of our heuristic on incompletely specified 
benchmarks. Table [3] reports a subset of our experiments. We have run the heuristic without the elimination 
of cubes covering only don't cares points from the solution under construction (DSOP-3 (a)), and with such 
elimination (DSOP-3 (b)). As the table clearly shows, the elimination of these cubes naturally produces 
better solutions in terms of size, but the computational time is much higher. 

In another series of experiments we compared our heuristic (with the third version of the optimization 
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Table 2: Comparison of five different variants of the DSOP minimization heuristic (SORT version: 
weight/dimension.) The size of the best DSOP representation computed for each benchmark is in bold- 
face. 

phase, and without elimination of cubes of don't cares only) with other DSOP minimization methods. We 
considered three techniques working, as ours, on explicit representation of cubes, and one method based on 
binary decision diagrams. The first algorithm [5 | sorts cubes in a minimal SOP according to their size, and 
compares the largest cube with all the others, starting from the smallest ones. In the next step, the second 
largest cube is selected and compared to all smaller ones, etc. As a last step, the cubes are merged wherever 
possible. The second algorithm, presented in lITSl . exploits the property of the most binate variable in a set 
of cubes to compute a DSOP form. The algorithm proposed in [2] enumerates all overlapping pairs of cubes 
in a SOP form, and builds a disjoint cover starting from the pairs of cubes with the highest degree of logic 
sharing. 

Finally, the third approach, presented in ||6], makes use of BDDs, exploiting the efficiency resulting from 
the implicit representation of the products. Observe in fact that a DSOP form can be extracted in a straight- 
forward way from a BDD, as different one-paths correspond to disjoint cubes. As the results presented in [6 1 
largely depend on the variable ordering of the underlying BDD, in |3 | an evolutionary algorithm has been 
proposed to find an optimized variable ordering for the BDD representation that guarantees more compact 
DSOP forms. 

Table HI reports a cost-oriented comparison among the different methods. The first three columns are as 
before. Columns four and five report the number of products in the PLA realization and in the SOP form 
heuristically minimized by ESPRESSO in the heuristic mode. The column labeled DSOP ESPR. shows the 
size of the DSOP computed running ESPRESSO with the option "-Ddisjoint" on the previously computed 
SOP form. The next five columns report the sizes, when available, of the DSOP forms computed with the 
methods discussed in Q, lfT5]| . lH, IS, and lH, respectively. Finally, the last column shows the size of the 
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Table 3: DSOP synthesis of incompletely specified benchmarks, without (DSOP-3 (a)) and with (DSOP-3 
(b)) elimination of cubes covering only don't cares. The size of the best DSOP is in boldface. 

DSOPs computed with our heuristic (third variant). 

As the table clearly shows, our method almost always generates smaller DSOP representations, and the 
gain in size can be quite striking, as for instance for the benchmarks alu4, clip and misexS. We have found 
only a few benchmarks where our approach compares unfavorably: 5xpl, cordic and inc. 

A time comparisons among all these different methods was not possible due to the partial absence of 
CPU times specification in the literature. 

6.2 Partial DSOP synthesis 

In order to test our partial DSOP synthesis algorithm, we have applied the heuristic to the classical ESPRESSO 
benchmark suite flSl with the following meaning. We have considered only benchmarks with don't cares, 
where the on set of the benchmark is the on set of sopD, and the don't care set of the benchmark is the don't 
care set of sopS. 

Table |5] reports a subset of our experimental results. The column labeled SOP shows the number of 
products in a SOP representation computed by ESPRESSO in the heuristic mode. The remaining three pairs 
of columns report the number of products and the corresponding synthesis time for the following three forms 
(all computed with the third version of the optimization phase, the dimension/weight sort version, and with 
the elimination of cubes covering don't cares only): 

1. DSOP: a DSOP for the original function, with the choice of don't cares performed by ESPRESSO in 
the heuristic mode. Each don't care point is covered at most once. 

2. P-DSOP (a): a partial DSOP for the original function, with the choice of don't cares performed by 
ESPRESSO in the heuristic mode. Don't care points are either eliminated or covered at least once. 

3. P-DSOP (b): a partial DSOP for the original function, where all the don't cares of the function are 
in play (they have all been covered during the first SOP minimization). Don't care points are either 
eliminated or covered at least once in the final form. 

Note that the results in the column SOP are better than ours because the resulting form is not disjoint. 

The table suggests that the best solution is the one relative to the choice of don't cares made by 
ESPRESSO. Moreover, it appears clearly from these results that the option of covering more than once 
the don't care points of the function (DSOP-3 (a)) gives better results, especially for big benchmarks. 
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16 


16 
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16 
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16 
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Table 4: Comparison with other techniques. The size of the best DSOP is in boldface. 

7 Conclusions and Future Work 

Deriving an optimal DSOP or partial DSOP representation of a Boolean function is a hard problem. This is 
why we have proposed a heuristic that has been implemented, tested, and compared with others. 

From the experimental results we conclude that exploiting SOP minimization for DSOP synthesis is 
a crucial idea. In fact, comparing our results with the ones in the literature we always obtain equal or 
smaller forms. We observe that the fact that SOP and DSOP problems are so close is not intuitive. In fact, 
we would have expected that efficient strategies to solve the two problems would be different since DSOP 
minimization appears to be much harder then SOP synthesis. Nevertheless, the experiments show that, 
starting from minimal or quasi-minimal SOP expressions, we can heuristically derive very compact DSOP 
forms. Moreover, from Table |5] we also infer that the choice of the don't cares, which are used as ones of 
the function, performed for the SOP minimization is nearly always the best choice also for DSOP synthesis. 
Therefore it would be interesting to further study the closeness of SOP and DSOP minimal forms both in 
theoretical and experimental way. 

It could also be worth studying the approximability of DSOP minimization with the aim of designing 
approximation algorithms instead of heuristics. In fact, while a /^-approximation algorithm yields a near- 
optimal solution, i.e. a solution whose cost C is < pC* where C* is the cost of an optimal solution [7], 
no prediction can be made on the result of a heuristic. Perhaps a first step in this direction would be 
understanding when our heuristic returns a DSOP whose cost is much higher then the cost of an optimal 
DSOP 
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Table 5: Partial DSOP synthesis, using the subset of don't cares selected by ESPRESSO-NON-EXACT (P- 
DSOP (a)) or all don't cares of the original function (P-DSOP (b)). 
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